Cosmography in F(G) modified gravity 

M. R. SetartS 
N. MohammadipouiH 

Department of Science, University of Kurdistan 
Sanandaj, IRAN. 



Abstract 

Investigating the accelerated expansion of the universe with cosmography 
is a best method to constraint cosmological models. In this work, in the 
F{G) modified gravity framework, we obtain equations of motion in a flat 
FRW metric. Then we reconstruct the present day values of F{G) and its 
derivatives with the cosmographic parameters on the only assumption that 
the universe is homogenous and isotropic on large scale. Also we 
investigate the conditions of cosmologically viable F(G) gravity models 
with the fiducial data set values. 
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1 Introduction 



Recent accelerated expansion of our universe is one the most significant 
cosmological discoveries over the last decade [H O [3lll] . In order to explain 
this interesting behavior, many theories have been proposed. Although it is 
widely accepted that the cause which drives the acceleration is the so called 
dark energy, its nature and cosmological origin still remain enigmatic at 
present. The simplest suggestion for dark energy is cosmological constant. 
But it suffers from two kind of problems [5] : fine tuning and coincidence 
problem. A dynamical scalar field with quintessence or phantom behavior 
is another proposal for dark energy ( for reviews see ). 

Apart from the mentioned possibility of the existence of dark energy, a 
prominent possibility is that the gravitational interaction is modified at (at 
least) cosmic scales. Amongst the latter, models generalizing the Einstein- 
Hilbert action have been proposed. The Einstein field equations of General 
Relativity (GR) were first derived from an action principle by Hilbert, by 
adopting a linear functional of the scalar curvature, R, in the gravitational 
Lagrangian density. However, there are no a priori reasons to restrict the 
gravitational Lagrangian to this form, and indeed several generalizations of 
the Einstein-Hilbert Lagrangian have been proposed. The simplest way for 
modification of GR is to replace Ricci scalar, R in Einstein- Hilbert action 
with a general function of the Ricci scalar which is well known as F{R) 
gravity [3 [8]. For this kind of modification, one assumes that the gravita- 
tional action may contain some additional terms which starts to grow with 
decreasing curvature and obtain a late time acceleration epoch. One of the 
other of modification of Einsteins gravity is Gauss-Bonnet (GB) modifica- 
tion. As a possibility the Einstein-Gauss-Bonnet gravity is low energy limit 
of the string theory is of particular interest because of its special features. 
The GB generalization adds quadratic terms, involving second order curva- 
ture invariants (specifically Gauss-Bonnet term is a topological invariant in 
four dimensions) to the Einstein-Hilbert Lagrangian [9j, (see also [TO]). As 
mentioned , in 4-dimension the GB term is trivial. Therefore, it is usually 
used from another form of modified gravity known as modified GB theory, 
e.g. see In this theory, modification is done by adding an arbitrary 

function F{G) into the Einstein-Hilbert Lagrangian. Also, in the present 
work, we focus on this model, but with an unknown function F{G) in hand 
that must be determined for expanding accelerated universe in spatially 
flat FRW geometry. It is worth noting that both dark energy models and 
modifled gravity theories have shown to be in agreement with the data. 

However a possible solution could be to come back to the cosmography 
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rather than finding out solutions of the Friedmann equations and testing 
them. In the present paper, we discuss the possibility to constrain F{G) 
gravity models using a cosmographic approach [12]. Cosmography relies 
on two crucial things: I) extracting the maximum amount of information 
from measured distances, like the luminosity distances of SNela, II) assum- 
ing that the universe is homogenous an isotropic on large scale. Recently, 
this approach was considered by using SN in pJJ, SN+GRBs in [l^ and 
SN+OHD+BAO in ^15j, where the current status of our universe can be 
read. 

This paper is outlined in the following manner: In Section[2]we briefly re- 
view modified Gauss-Bonnet gravity, for self-completeness and self-consistency, 
and present the respective field equations. In section[3]we introduce the basic 
notions of the cosmographic parameters then, write down all the quantities 
relevant of the F{G) models (are presented in previous section) in terms of 
the cosmographic parameters. Inserting cosmographic parameters of ACDM 
model into the expressions of the F{G) quantities, we will investigate in sec- 
tion m Section [5] contains the main result of the paper demonstrating how 
the F{G) derivatives can be related to the cosmographic parameters and we 
use these results and constraints on the cosmographic parameters to derive 
the cosmologically viable of F{G) models and model independent estimates 
of the present day values of the F{G) derivatives. Finally, in Section [6l we 
present our conclusions. 



2 Field equations and F{G) modified gravity 

We start with the 4-dimensional action in F{G) gravity 

dx^ + S^, (1) 

where IZ is Ricci scaler curvature, F{G) is a generic function that must be 
determinecH, k"^ = Svr^ and Sm is the matter action. 

Varying the action with respect to the metric, one can obtain the field 
equation as 

G^u + s\TZ^piia + TZpi/Qo-fi + — {gfii/go-p ~ gfiagi/p)TZ ~ Ttpcjgup, (2) 

-n^ugap + n^agup]VPV''F'{G) + (GF'(G) - F)gp, = , 
^Indeed, G is topological invariant G = TZ^ - ATZcfiTl"^ + TZafj-ysTZ"'^^^ . 



S 



1 

fc2 



7^ 



+ F(G) 
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where G^^ the Einstein's tensor, is the energy-momentum tensor of 
the matter, k"^ = Svr^ = 1 and prime denotes the ordinary derivative with 
respect to G. For spatiahy flat Robertson- Walker metric 

ds^ = -dt^ + a(t)2dx2, (3) 

we have 

n = 6{H + 2H^), G = 24H^{H + H^) , (4) 

where H is the Hubble parameter and dots denotes the time derivative. 
Considering that the universe is filled by pressureless dust matter, under 
the assumption of a flat universe the Eq.Q read 

3H^ = GF'{G)-F-2AH^F'{G) + p„„ (5) 

-2H = -SH'^F'iG) + WHHF'{G) + 8H^F'{G) + pm, (6) 
and also continuity equations for the densities pm satisfy the following form 

Pm + 3Hpra = 0. (7) 

Considering F{G) gravity as a factor that, is causing the accelerated expan- 
sion of the universe, we can rewrite the energy density and pressure of this 
dark energy as follows 

PG = GF\G) -F- 2AH^GF"{G), (8) 



PG = [{IQH^G + 16HHG + 8H^G)F"{G) (9) 
+8H^G'^F"'{G) - GF'{G) + F]. 

The continuity equation for pG leads to 

PG + 3HipG + Pg) = 0; (10) 

with the equation of state as 

im'^G + IQHHG - SH^G)F" + SH^G^F'" , , 

= — 1 H : . (11) 

GF' -F - 2AH3GF" 

Here, quantities pG and log depends on H, G, F{G) and there derivatives. 
These quantities could be related to the observed acceleration of the uni- 
verse. 
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3 F(G) gravity vs cosmography 



In this section, firstly defined cosmographic parameters then, instead of solv- 
ing Eq.(l6]) for given values of the boundary conditions, we obtain the relation 
between the present day values of F{G), its derivatives and cosmographic 
parameters. 

A. Cosmography parameters 

The cosmographic parameters that are proportional to the coefficients 
Taylor series expansion of the scale factor with respect to the cosmic time 
defined as 

H(t) = (12) 

a at 

11 cPa , , 

= -IPaW 

11 d'^CL , , 

11 / N 

= /Fad?' <^'=> 

which are usually referred to as the Hubble, deceleration, jerk, snap and lerk 
parameters, respectively. For example, the present value of Hubble param- 
eter describes the present expansion rate of our universe, and a negative 
value of go means that our universe is undergoing an accelerated expansion 
while, jo allows to discriminate among different accelerating models. We can 
rewrite the Hubble parameter and its derivatives in terms of cosmographic 
parameters as 

i? = -F2(l + g), (17) 
if = /73(j- + 3g + 2), (18) 



//^s - 4j - 3g(g + 4) - 6, (19) 



H'^il - 5s + 10{q + 2)j + 30(g + 2)q + 24], (20) 
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where a dot denotes derivative with respect to the cosmic time. Using 
Eqs. (|17p -()20 p we can rewrite the Gauss-Bonnet term and its derivatives in 
terms of cosmographic parameters. 

B.F(G)-derivatives and cosmography 

Rather than solving Eq.([6]) and choosing F{G) by given values of bound- 
ary conditions, we reconstruct the present day value of F{G) and its deriva- 
tives with the present time of the cosmographic parameters {Hq, qo,jo,, Sq, Iq)- 

Considering Eq.(j3|) and differentiating with respect to t, one can obtain 
the following equations 

G = 24[2(2F3 + HH)H + H^H], (21) 



G = 24[2(6i?2 + H)H^ + 2{2H^ + 3HH)H + H^-^r], (22) 



cfiG 

24[24:HH^ + 6{6H'^H + 2H'^ + HH)H (23) 



dt2, 



Since we will reconstruct the present day value of F(G) and its derivatives, 
then rewrite Eqs.(l5]), ([6]) at the present time t = tg as 



2 _ GoF'jGo) F{Go) 3 ■ „ 2o 

-"0 — ^ ^ oHqUo^ [Uq) + HQ\lom, 



(24) 



Ho = i^HlGo - 8H0H0G0 - 4HlGo)F"{Go) (25) 

-4i/gGo F"'{Gq) ^Oom, 

where pm = Pmi't = to)o~^ = S-fTgriomO"^! with fiom the present day matter 
density which is obtained from the matter continuity equation. 

We will assume that F{G) may be well approximated by its third order 
Taylor expansion in (G — Go) following as 

F{G) = F{Go) + F'{Go){G - Go) + ^F"(Go)(G - Gof (26) 
+JF'"(Go)(G-Go)^ 
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with Go being the current value of G. Here we assume ' = for n > 4. 
In order to get the third derivative of F{G), we differentiate of Eq.([6]) with 
respect to t and evaluate in (t = to)- One finally obtain 

Ho = ^Qora + mHoHo - 2Ho - 2HoHo)Go + (27) 

{Hq — 4HoHo)Go — Hq—-j^]F"{Go) — 

4[3HlGoGo - {Hi - 4HoHo)Go^]F"'{Go). 

Evaluate the present day value of {H, H, ^) Eqs.(II7])-(I5D]) and sub- 
stitute these equations into Eq.([3]) and Eqs. ([2T]) - (p3]) to obtain G and its 
derivatives at the present time: 

Go = -24^0^90, (28) 
Go = 24Fo^(2gg + 3(?o+io), (29) 
Go = 24Hl[iqo + l)(-2g2 - 12go - 6 jo) - + so], (30) 

= 24Hl^[{qo + l)(36g2 + 60go - 8so) + 12gg (31) 

+/o -S0 + 6io(2Qo + llgo + JO + 6)]. 

Using Eqs.([17|-(I20]) and by inserting Eqs.(l28l)-(l3l|) into Eq.([25]) and Eq.(l27D 
finally, we reconstruct the present day value of the second and the third 
order derivatives of F{G) with respect to cosmographic parameters at z = 
as 

F"{Go) _ Ao2 + Bo2^0m /„„x 

(24i76)-i - Go ' ^^'^ 



F'"(Go) Ao3 + Bo^a 



Om 



(24i/5)- 



Gn 



(33) 



with 



Ao2=4{2ql+3qo+jo)ijo 



21 jo - 39goio 



20q^oh 



(34) 



-75q^ + 3qoSo 



14'Zo' 



105gg + 3so-45go), 
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Bo2 = -6{2ql + 3qo + io)(16jo + Uqoh + 29gg - 3so - 2ql + 30go), (35) 



Ao3 = 48so - 228jo - - 420(7o - 888(7^ - 408goJo - 762q^ - 298g^io (36) 



+4j2 - 404q^ - lUq^jo - 48g^ + 48q^So - 4joSo + 60goSo - ^qok 



Bo3 = 288J0 + 6/0 - 66so + 50Aqo + I242ql + 690(7oio (37) 
+720gg + 252g2 jo + ASjl - 72qoSo + 72ql 



Co = -16(2g^ + 3qo + jo){8j^ - 132j^ - 219j^qo - llSq^^j^ (38) 

+36joSo + - 570gojo - 148goio - 1239goJo - Ssq 
-852g3jo + 34jogoSo - lQ26ql - GSOq^ + 84^2^0 + 3^70/0 
+72^oSo - 1419(7^ - 60ql + 2qllo - 468g^ + Sq^so). 



Insert Eg . (1321) into Eg . (1241) and we get to relation between the present day 
value of F{G) and its first derivative 

F'{Go) , F(Go) _ ^ (Co + 8aoAo2) + {8aoBo2 - Co)Oom 



(24i/2g4)-l ^2 

where 

ao = (2g2+3go+jo)- (40) 

Egs. (f32]l - (fl0]l make it possible to estimate the present day values of the 
second, the third derivatives of F(G) and the relation between F(G) and its 
first derivative as function of the Hubble constant Ho and the cosmographic 
parameters {qo, jo, sqJq) provided a value for the matter density parameter 
Horn is given. Although, the value of Clom is usually the result of fitting 
a given dataset in the framework of an assumed dark energy scenario, but 
different models all converge towards the concordance value ^lom ~ 0.25 
which is also in agreement with astrophysical (model independent) estimates 
from the gas mass fraction in galaxy clusters. Indeed, Egs. ()32p - ()40p are ideal 
tools to testing F(G)-gravity models with observational dataset . 

4 The ACDM model 

We can test the reliability by comparing result are obtained in previous 
section with spatially flat ACDM model. The cosmographic parameters for 
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the ACDM model with equation of state parameter oj = —1 at the present 
time read 



qo = - 




(41) 



jo = l 



(42) 



(43) 



Iq — 1 + Sfiom + -Z-^Om- 



(44) 



Inserting Eqs.(|ll])-(I111) into Eqs.(l32]) and ^ gives 



F"(Go) = F'"(Go) = 0. 



(45) 



From Eq. (j3ip and the previous equations in Eqs. (|4ip - ()44p one can obtain 



The equation of state for these models defined by Eq. (fTT]) that, by us- 
ing Eq. ([l5]) reduce to ujg = —1- Thus, we obviously conclude that the 
R + F{G) theory having exactly the same cosmographic parameters as the 
ACDM model with equation of state ujg = -1 if GoF'(Go) - F{Go) = A. 
If we assume GoF'{Go) = 2A then F{Go) = A and the action Eq.([T]) reduce 
to GR + A + matter. 



5 Constraining F{G) parameters 

In order to constrain the model by cosmography, i.e. to estimate the function 
GoF'(Go) - F{Go), F"{Go) and F"'{Go) at the present time, we need to 
constrain observationally the cosmographic parameters by using appropriate 
distance indicators. Moreover, we must take care that the expansion of the 
distance related quantities in terms of {qo,jo, sq,Io) closely follows the exact 
expressions over the range probed by the data used. 



F{Go) 



A + GoF'(Go). 



(46) 



A. Observational constraints 
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X 


XBF 


(x) 




68% CL 


95% CL 


h 


0.718 


0.706 


0.706 


(0.693, 0.719) 


(0.679, 0.731) 


qo 


-0.64 


-0.44 


-0.43 


(-0.60, -0.30) 


(-0.71, -0.26) 


jo 


1.02 


-0.04 


-0.15 


(-0.88, -0.90) 


(-1.07, 1.40) 


So 


-0.39 


0.18 


0.02 


(-0.57, 1.07) 


(-1.04, 1.78) 


k 


4.05 


4.64 


4.54 


(2.99, 6.48) 


(1.78, 8.69) 



Table 1: Constraints on the cosmographic parameters by jointly fitting the 
Union2 SNela sample and the BAO data. Columns are as follows : 1. pa- 
rameter id; 2. best fit; 3., 4. mean and median from the marginalized 
likelihood; 5., 6. 68 and 95% confidence ranges. 

Taking SNela and a fiducial ACDM model as a test case, one has to check 
that the approximated luminosity distanc^ deviates from the ACDM one 
less than the measurement uncertainties up to z ~ 1.5 to avoid introducing 
any systematic bias. Since we are interested in constraining ((?0) JO; ^o)) 
we will expand the luminosity distance up to the fifth order in z which 
indeed allows us to track the ACDM expression with an error less than 
1% over the full redshift range. We have checked that this is the case also 
for the angular diameter distance Da = Di{z)/{1 + z)"^ and the Hubble 
parameter H{z) which, however, we expand only up to the fourth order to 
avoid introducing a further cosmographic parameter. In order to constrain 
the parameters {h,qo,jo,so,lo), we consider the dataset that is used in [12] 
and are summarized in Table[TJ 

In order to translate our constraints on the cosmographic parameters on 
similar constraints on F{G) and its derivatives, we should just use Eqs. (|32p . 
(|33p and (I39p evaluating them along the final coadded and thinned chain 
and then looking at the corresponding histograms. We define quantities for 
shortness as 

/o = GoF'(Go)-F(Go), /a = F"(Go), f, = F"'{Go). (47) 
*See [12] for the analytical expression. 
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X 


XBF 


(x) 


^med 


68% CL 


95% CL 


/o 


1.139 


0.921 


0.905 


(0.859, 0.854) 


(0.765, 1.801) 


/2 


0.0027 


0.0050 


0.0050 


(0.0030, 0.0045) 


(0.0038, 0.0321) 


/3 


-0.0231 


-0.0020 


-0.0019 


(-0.0011, -0.0010) 


(-0.0015, 0.0734) 




-1.0238 


-0.8542 


-0.8454 


(-1.0139, -0.0718) 


(-1.1334, -0.9849 ) 



Table 2: Constraints on quantities /j and equation of state values from the 
Markov Chain for the cosmographic parameters. Columns order is the same 
as in Table[TJ 

We constrain these quantities by setting the value of f^om along the 
chain using r^om = 0Jmh~'^ with the physical matter density oom = 0.1329 in 
agreement with the WMAP7 data. For each value of h along the chain, we 
fix = ^ml^ having neglected the error on uj^ since it is subdominant 
with respect to the one on h. We stress that, although the fiducial value for 

has been obtained for a ACDM model, it should be unchanged for any 
model which reduces to the GR + matter domination at the CMBR epoch 
as is our case. 

We finally analyize the present day value of quantities summarized in 
Table[2j It is worth noting that the /j values become smaller as the order 
n of the derivative increases and it is in agreement with our assumption in 
Eq. (l26p . Note also that, mean and median from the marginalized likelihood 
values of {f2, fs) are almost equal and are not quite different from their 
median values. As a further remark, we note that the confidence ranges 
of constraints on fo, f2, fs (in particular, for 95%) become larger as the 
order i of the derivative increases. Which caused the degeneracies among 
(90)ioi So, /q) and the nonlinear behavior of the relations /«. 

B. Observational constraints on cosmologically viable F{G) 

gravity models 

In this section, we translate the constraints on the cosmographic param- 
eters on similar constraints on the present day values of conditions cosmo- 
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X XBF {x) Xmed 



68% CL 



95% CL 



/2 0.0027 0.0050 0.0050 (0.0030, -0.0050) (0.0038, 0.0321) 
H^f2 0.0004 0.0006 0.0006 (0.0003, -0.0007) (0.0004, 0.005) 



Table 3: Constraints on Mq values from the Markov Chain for the cosmo- 
graphic parameters. Columns order is the same as in Tabled) 



logically viable F{G) models. De Felice et al. [16| have investigated the 

viability of cosmological F(G) gravity models. They started work with our 
action Eq.([T]) in a spatially FRW background with the metric signature (- 
, +, +, +) and showed that the viable of this models need to satisfy the 
following conditions: 

• (1) F{G) and its derivatives F'{G), F"{G),... are regular. 

• (2) F"{G) > for G < Gi and F"{G) approaches +0 in the limit 
|G| ^ oo. 

• (3) < HfF"{Gi) < 1/384 ~ 0.0026 is the stable condition at the de 
Sitter point. 

• (4) < HfF"{Gi) < 1/600 ~ 0.0017 corresponds to a stable spiral 
(damping with oscillations). 

Here index 1, in this context, denotes the calculated values of this quan- 
tities in de-Sitter point. In order to investigate the above conditions with 
the cosmographic parameters at the present time, we use Table [T] and eval- 
uate Mq = HqF" {Gq). Finally, the present day values of /2 = F"{Go) and 
Mq summarize in Tabl^ 

The observational constraints summarized in Tabl^ shows that, the 
present day value of /2 = F"{Gq) > and it is in agreement with the 
similar results in [16] to have a stable de - Sitter point in F{G) models. We 
also note that the values Mq = HqF"{Gq) in Tabl^ corresponds to a stable 
spiral de - Sitter point (damping with oscillations). 




C. Cosmography VS F{G) models 
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Up to now, for a generic function of F{G), we obtain observational con- 
straints on F{G) and its derivatives from cosmographic parameters. Now, 
we check the viabihty of a given F{G) model with constraints are provided 
in previous section without the need of explicitly solving the field equations 
and fitting data. 

As an example, we assume the following model as|17] 

= aG" + ^G In G, (48) 



where a, /3 and n are constants. Using quantities Eq. ()17|) and Eqs. (f32]) . p3]) 
and ([5U|1 . one can obtain the following expressions 

/o = a(n-l)G^ + /3Go, (49) 
/2 = an(n - 1)G^2 + /3Go \ (50) 
h = an{n-l){n-2)G^-'-PG,\ (51) 

Solving the first and second equations with respect to (a, f3) leads to 



/o ~ 
(n-1) 



« = - L T2 go"^ (52) 



/3 = ^^0-^0/2 ^-1^ (53) 
n — 1 

Inserting a and /? into the third equation of Eq.(51), we obtain two set 
solutions for n. 



-P(/0i f2, fs, Gq) ± Q(/o, /2, /s, Gq) , , 

(Jo, h, (j^oj 



with 



P(.foj2j3,Go) = /3 + 2G0V2-G0V0, (55) 
r(/o,/2,Go) = /2Go^-/oG'o^ (56) 

Q(/o, /2, /3, Go) = Go Vo + 4Go VI + fi - 4/0/2G0 ^ (57) 
+2/0/3G03-4/2/3G0I. 

We can see that values of (a,/3,n) depends on the values of {qo, jo, sq^Iq) 
then, one may obtain one, two or any acceptable solution, i.e. real positive 
values of n. If the final values of (a, /?, n) are physically viable, we can 
conclude that the model in Eq. (j48p is in agreement with the data giving the 
same cosmographic parameters inferred from the data themselves. 
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6 Conclusion 



Cosmography is a ideal way to give a picture of the observed universe con- 
sidering only assumptions that, the universe is homogenous and isotropic on 
large scale. In the present paper we have discussed the cosmography of a 
modified GB model with an unknown function of topological invariant GB. 
At first we have defined cosmographic parameters, then instead of solving 
equations of motion for given values of the boundary conditions, we have 
obtained the relation between the present day values of F{G), its deriva- 
tives and cosmographic parameters. We have assumed that F{G) may be 
well approximated by its third order Taylor expansion in (G — Go). Then we 
have reconstructed the present day value of the second and the third order 
derivatives of F{G) with respect to cosmographic parameters at z = by 
Eqs. (f32]) . ([33]) . We also found a relation between the present day value of 
F{G) and its first derivative. This recent relation and also the second and 
the third order derivatives of F{G) are functions of the Hubble constant Hq 
and the cosmographic parameters {qo,jo,so,lo), and $7om • 
In contrast to the cosmography investigations in f{R) and /(T) [12] here 
we note that, in our model do not have to considering this assumption 
that Geff{z = 0) = Gn — )• /'(O) = 1. We transfer the constraints on 
(-ffo; 90; jO) •50, /q) and Qom into similar ones for fi quantities and summa- 
rized them in Table[2l It is worth noting that, these constraints (comes out 
from fi) showed that, the validity of our assumption in the Taylor expan- 
sion of F{G) which is = for n > 4. On the other hand, it is clear 

that the increasing order of expansion shifts away from the ACDM fiducial 
values. Then, it is possible recover the ACDM with the above assumption 
in F{G) gravity. 

As a further remark, we have investigated the conditions of cosmologi- 
cally viable F{G) gravity models with the fiducial values in TableO These 
values were in agreement with the viable conditions of this model that the 
authors of [16] investigated, so that the stable condition is /2 = F"{G) > 
and < HfF"{Gi) < 1/600 ~ 0.0017 corresponds to a stable spiral (damp- 
ing with oscillations). 
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